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Abstract 

The angular bispectrum of spherical random fields has recently gained 
an enormous importance, especially in connection with statistical infer- 
ence on cosmological data. In this paper, we provide expressions for its 
moments of arbitrary order and we use these results to establish a mul- 
tivariate central limit theorem and higher order approximations. The re- 
sults rely upon combinatorial methods from graph theory and a detailed 
investigation for the asymptotic behaviour of Clebsch-Gordan coefficients; 
the latter are widely used in representation theory and quantum theory 
of angular momentum. 
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1. INTRODUCTION 



Let T(9, (f) be a random field indexed by the unit sphere S 2 , i.e. < 9 < tt 
and < (p < 2ir. We assume that T(9, ip) has zero mean, finite variance and 
it is mean square continuous and isotropic, i.e. its covariance is invariant with 
respect to the group of rotations. For isotropic fields, the following spectral 
representation holds in mean square sense (Yaglom (1986), Leonenko (1999)): 



oo 



T(B,<p)=2^ ^ a lm Y lm (9,<p) . (1) 

£=1 m= — l 



*I am very much grateful to M.W.Baldoni for many discussions and explanations on the 
role of Clebsch-Gordan coefficients in representation theory. Usual disclaimers apply. 
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Here, we have introduced the spherical harmonics (see Varshalovich, Moskalev 
and Khersonskii (hereafter VMK) (1988), chapter 5), defined by 



Yi m (9,<p) 




Pim (cos 9) exp(jTO<^) , for m > 



Y lm (9,<p) 



(-ir^-m(^) ; form<0 



where the asterisk denotes complex conjugation and Pi m (cos9) denotes the as- 
sociated Legendre polynomial of degree l,m, i.e. 



The triangular array {ai m } represents a set of random coefficients, which can 
be obtained from T(0, ip) through the inversion formula 



see for instance Kim and Koo (2002), Kim, Koo and Park (2004) for a review 
of Fourier analysis on S 2 . The coefficients ai m are complex- valued, zero-mean 
and uncorrelated; hence, if T(9, ip) is Gaussian, they have a complex Gaussian 
distribution, and they are independent over I and m > (although a; ; _ m = 
(— l)" l af m ), with variance E\ai m \ 2 = C\, m = 0, ±1, ±Z. The index I is usually 
labeled a multipole; approximately, each multipole corresponds to an angular 
resolution given by n/l. 

The sequence {C{\ denotes the angular power spectrum: we shall always 
assume that Ci is strictly positive, for all values of I. As well-known, if the field 
is Gaussian the angular power spectrum completely identifies its dependence 
structure. For non-Gaussian fields, the dependence structure becomes much 
richer, and higher order moments of the aim's are of interest; this leads to the 
analysis of so-called higher order angular power spectra. 

The analysis of spherical random fields has recently gained momentum, due 
to strong empirical motivations arising especially (but not exclusively) from 
cosmology and astrophysics. In particular, an enormous attention has been 
drawn by issues connected with the statistical analysis of Cosmic Microwave 
Background radiation (CMB). CMB can be viewed as a snapshot of the Uni- 
verse approximately 3 x 10 5 years after the Big Bang (Peebles (1993), Peacock 
(1999)). A number of experiments are aimed at measuring this radiation: we 
mention in particular two satellite missions, namely WMAP by NASA, which 
released the first full-sky of CMB fluctuations in February 2003, with much 
more detailed data to come in the years to come, and Planck by ESA, which 
is due to be launched in Spring 2007 and expected to provide maps with much 
greater resolution. Over the next ten years, an immense amount of cosmolog- 
ical information is expected from these huge data sets; at the same time, the 



m = 0,1,2, , I = 1,2,3, 



aim = 




T(6,<p)Y lm {$,¥>) sin 6d6dip , m = 0, ±1, ±1 , I = 1,2, 



(2) 
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analysis of such data sets posits a remarkable challenge to statistical methodol- 
ogy. In particular, several papers have focussed on testing for non-Gaussianity 
by a variety of nonparametric methods (to mention a few, Dore, Colombi and 
Bouchet (2003), Hansen, Marinucci and Vittorio (2003), Park (2004), Marin- 
ucci and Piccioni (2004), Jin et al. (2005)). The majority of efforts has focussed 
on the angular bispectrum, which is considered an optimal statistic to verify 
the accuracy of the so-called inflationary scenario, the leading paradigm for 
the dynamics of the Big Bang. See for instance Phillips and Kogut (2000), 
Komatsu and Spergel (2001), Bartolo, Matarrese and Riotto (2002), Komatsu 
et al. (2002,2003), Babich (2005) and several others; a review is in Marinucci 
(2004). 

The angular bispectrum can be viewed as the harmonic transform of the 
three-point angular correlation function, whereas the angular power spectrum 
is the Lcgcndre transform of the (two-point) angular correlation function. Write 
fli = (9i, tp t ), for i = 1,2, 3; we have 

oc 

ST(n x )T(n 2 )T(0 3 ) = E B hhh m3Y ^ i^) Y hm 2 (n 2 )Y l3m3 (fi a ) 

li,l 2 ,l3 = l m 1 ,m2,m 3 

(3) 

where the bispectrum -B™^ 2 ™ 3 is given by 

llhh ~ £j \ a hm 1 0,l 2 m 2 a l3'm 3 ) ■ l*J 

Here, and in the sequel, the sums over m% run from —U to U, unless otherwise 
indicated. Both J3J and (@J are clearly equal to zero for zero-mean Gaussian 
fields. Moreover, the assumption that the CMB random field is statistically 
isotropic entails that the right- and left-hand sides of © should be left unaltered 
by a rotation of the coordinate system. Therefore S™^™ 2 ™ 3 must take values 
ensuring that the three-point correlation function on the left-hand side of J2J 
remains unchanged if the three directions fli , ft 2 and il 3 are rotated by the same 
angle. Careful choices of the orientations entail that the angular bispectrum of 
an isotropic field can be non-zero only if li < lj + Ik for all choices of i, j, k = 
1, 2, 3; li + h + h is even; and mi+m2 + ni3 = 0. More generally, Hu (2001) shows 
that a necessary and sufficient condition for i?™^™ 2 ™ 3 to represent the angular 
bispectrum of an isotropic random field is that there exist a real symmetric 
function of l\, h, hi which we denote bi^i^, such that we have the identity 

R ro ,T 2m3 - GTfT m3 Ki2h ; (5) 

bi t i 2 i 3 is labeled the reduced bispectrum. In (JSJ we are using the Gaunt integral 
^ 2m3 , defined by 

7T r2ir 



GZT m3: = / Y limi (9,<p)Y l2m2 (e,ip)Y hms (6,<p)sin6d<pd6 



(2/ 1 + l)(2/ 2 + l)(2; 3 + l)\ 1/ V h h k \ (h h I. 



47r J \ J \ mi m 2 m 3 



'3 
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where the so-called "Wigner's 3j symbols" appearing on the second line are 
defined by (VMK, expression 8.2.1.5) 



h h h 
mi m 2 m 3 



(-1) 



(h + l2-h)Kh-l2 + h)Kh-h + hy. 



1/2 



(h + h + h + iy. 

(l 3 +m 3 )\(l 3 -m 3 )\ l 1 / 2 



(h + mi y.{h - mi y.(i 2 + m 2 y.{i 2 



m 2 ) 



E 



{-l) z {h + h + mx- z)\(h ~ mi + z)\ 
z\(h +h-h- z)\{h + m 3 - z)\(h -h-m 3 + z)\ 



where the summation runs over all z's such that the factorials are non-negative. 
Note that the Wigner's 3j are invariant with respect to permutations of the 
pairs (Zj,mi). 

In view of (J5J, the dependence of the bispectrum ordinates on mi,m 2 ,mz 
does not carry any physical information if the field is isotropic; hence it can be 
eliminated by focussing on the angular averaged bispectrum, defined by 



h h h 
mi m 2 m 3 



h h h 

Bhi 2 i 3 ■■= E E 

mi- — li m.2 — — 12 mz — — lc 

(2^ + l)(2^ + l)(2l 3 + l) \ 1/2 / h i 2 i 3 
Air y ^000 

where we have used the orthogonality condition 

£1 h I3 



(G) 



E E E 

mi= — h m 2 = — (2 ni!=- (3 



h h h 
mi m 2 m 3 



= 1 



The minimum mean square error estimator of the bispectrum is provided by 
(Hu 2001) 



E E E 

mi= — ii 7712— — 7 2 7773— —£3 



h h h 
mi m 2 m 3 



(^ 7 l 777l ^Z 2 7772 ^^37773 ) ■ 



The statistic Bi t i 2 i 3 is called the (sample) angle averaged bispectrum; for any 
realization of the random field T, it is a real-valued scalar, which does not 
depend on the choice of the coordinate axes and it is invariant with respect to 
permutation of its arguments Zi, l 2 , l 3 . 

Under Gaussianity, the bispectrum can be easily made model- independent, 
namely we can focus on the normalized bispectrum, which we define by 



(-1) 



(h+h+h)/2_ 



B, 



lll'2h 



(7) 



The factor (— i)('i+'2+'3)/ 2 j g usua vjy no t included in the definition of the normal- 
ized bispectrum; it corresponds, however, to the sign of the Wigner's coefficients 
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for mi = ni2 = 7713 = 0, and thus it seems natural to include it to ensure that 
Il t l 2 l 3 and 6; 1 ; 2 ; 3 share the same parity (see ||BJ|). 

In practice Ii t i 2 i 3 is unfeasible because C\ is unknown. A natural estimator 
for Ci is 



which is clearly unbiased (see also Arjunwadkar et al. (2004)). Thus I^i^ can 
be replaced by the feasible statistic 



Although the bispectrum has been the object of an enormous attention in 
the cosmological literature, very few analytic results are so far available on its 
probabilistic properties. In a previous paper (Marinucci (2005)), we established 
some bounds on the behaviour of its first eight moments, and we used these 
results to establish the asymptotic behaviour of some functionals of the bispec- 
trum array; such functionals were proposed to build nonparametric tests for 
non-Gaussianity. The results on these moments where established by means 
of a direct analysis of cross-products of Wigner's 3j coefficients; this analysis 
was performed by means of explicit summation formulae originated from the 
quantum angular momentum literature (see VMK for a very detailed collection 
of results). In the present paper, these results are made sharper and extended 
to moments of arbitrary orders by means of a more general argument. More 
precisely, we show how it is possible to associate to higher moments the co- 
efficients of unitary matrices transforming alternative bases of tensor product 
spaces generated by the spherical harmonics. We are then able to obtain combi- 
natorial expressions for (cross-)moments of arbitrary orders. These results are 
then exploited to obtain multivariate central limit theorems and higher order 
approximations. It should be noted that the asymptotic theory presented in 
this work is of a fixed-domain type, a framework which has become increasingly 
popular in recent years, see for instance Stein (1999) or Loh (2005). 

The structure of the paper is as follows: in Section 2 we review some basic 
combinatorial material on diagrams and graphs; in Section 3 we present our 
general results on moments and we exploit it to obtain multivariate central limit 
theorems for the angular bispectrum with known or unknown C;; in Section 4, 
we discuss higher-order approximations. 



We shall review here some elementary notions from graph theory, which 
is widely used in physics when handling Wigner's 3j coefficients (see VMK, 






2. DIAGRAMS AND GRAPHS 
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chapter 11). Take i = l,...,p and j = 1, q t , and consider the set of indexes: 



we stress that the number of columns qi need not be the same for each row i. A 
diagram 7 is any partition of the elements of T into pairs like {(ii, ji), («2, J2)} : 
these pairs are called the edges of the diagram. We label T(T) the family of these 
diagrams. We also note that if we identify each row ik with a vertex (or node), 
and view these vertexes as linked together by the edges {(ik,jk), {h',jk')} = 
ikik', then it is possible to associate to each diagram a graph. As it is well 
known, a graph is an ordered pair (I, E) where I is non-empty (in our case the 
set of the rows of the diagram), and E is a set of unordered pairs of vertexes 
(in our cases, the pair of rows that are linked in a diagram). We consider only 
graphs which are not directed, that is, (M2) and (1211) identify the same edge; 
however, we do allow for repetitions of edges (two rows may be linked twice), 
in which case the term multigraph is more appropriate. A graph carries less 
information than a diagram (the information on the "columns" , i.e. the second 
element jk, is neglected) but it is much easier to represent pictorially. We shall 
use some result on graphs below; with a slight abuse of notation, we denote the 
graph 7 with the same letter as the corresponding diagram. 
We say that 

a) A diagram has a flat edge if there is at least a pair {(ii, Ji), (*2, J2)} such 
that i\ =12', we write 7 <E Tp(T) for a diagram with at least a flat edge, and 
7 G ^-p{T) otherwise. A graph corresponding to a diagram with a flat edge 
includes an edge ikik which arrives in the same vertex where it started; for 
these circumstances the term pseudograph is preferred by some authors (e.g. 
Foulds (1992)). 

b) A diagram 7 € Y-p(T) is connected if it is not possible to partition the 
fs into two sets A, B such that there are no edges with i\ € A and *2 € B. We 
write 7 G Tc{T) for connected diagrams, 7 e T^(T) otherwise. Obviously a 
diagram is connected if and only if the corresponding graph is connected, in the 
standard sense. 

c) A diagram 7 e IV(T) is paired if, considering any two set of edges 
{(n, ji), (t2,h)} and {(« 3 ,.73), {U, ji)} , then i\ = i 3 implies i 2 = ir, in words, 
the rows are completely coupled two by two. We write 7 G Tp(T) for paired 
diagrams. 




G 



Figure I: 7 G Tp(T) 



d) We shall say a diagram has a k-loop if there exist a sequence of k edges 
{(ii.ii), («2, J2)} , {(ifc, jk), (ifc+i, ife+i)} = («i«2), -, (ifcifc+i) 

such that zi = ifc+i; we write 7 € T^^^T) for diagrams with a fc-loop and no 
loop of order smaller than k. 

Note that Tp(T) = T L ^(T) (a flat edge is a 1-loop); also, we write 



for connected diagrams with k-loops, and r cjjj^{T) for connected diagrams 
with no loops of order k or smaller. For instance, a connected diagram belongs 
to r c -^y(T) if there are neither flat edges nor two edges {(«i, ji), {12,32)} and 
{(*3i J3), (*4)J4)} such that i\ — 13 and 12 = ii] in words, there are no pairs of 
rows which are connected twice. 

e) A tree is a graph with no loops (written 7 £ Tr(T)). 

Graphs and diagrams play a key role to evaluate the behaviour of moments 
of the bispectrum; to this issue we devote the next section. 



3. A CENTRAL LIMIT THEOREM FOR THE BISPECTRUM 

In this section, we shall investigate the behaviour of the higher order mo- 
ments for the normalized bispectrum Q, under the assumption of Gaussianity; 
to this aim, we define 



rcL(k)(T) = r c (T)nr L{k) (T) 




1 for h <l 2 < h 
2 for h = I2 < h or h < h = h 
6 for l\ = I2 = h 



here and in the sequel, S b a denotes Kronecker's delta, that is S b a = 1 for a = b, 
zero otherwise. 
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Under Gaussianity, it is obvious that the expectation of all odd powers of 
Ihhh is zero. To analyze the behaviour of even powers, we first recall that, 
for a multivariate Gaussian vector (xi, ...,X2k), we have the following diagram 
formula 

E(xi x x 2 x ... x x 2 k) = /XEx^Xjz) x ... x {Ex i2k „ 1 x l2k ) , (9) 

where the sum is over all the (2k)\/{k\2 k ) different ways of grouping [x\, ...,X2k) 
into pairs (see for instance Giraitis and Surgailis (1987)). Even powers of 7; 1 j 2 ; 3 
yield even powers of the a/ m 's, which have a complex Gaussian distribution, 
weighted by Wigner's 3j coefficients. 

In the sequel, unless otherwise specified, we rearrange terms so that In < 
l%2 < l%3 for all i. 

Theorem 3.1 Assume that (^1,^2,^3) 7^ {k'l, h'2, h'3) whenever i ^ i'. There 
exist an absolute constant K P1 ,_ PI such that, for pi > 1 , i = 1, 

U=l ) i=l 

always, where (2p — 1)!! = (2p — 1) X (2p — 3) x ... x 1. 

Remark 3.1 The condition that (^1,^2,^3) ^ {k<iji'2,h'3) whenever % ^ i' is 
merely notational and entails no loss of generality; indeed, whenever (Z<i, Zj2, ^3) = 
(h'i, h'2, h'3) it suffices to identify the two indexes and change the values of pi 
accordingly. 

Proof The proof is lengthy and computationally burdensome, so before we 
proceed we find it useful to sketch heuristically its main features. The first 
step is to notice that, in view of 10, higher order moments can be associated 
with sums over all possible graphs configurations of cross-products of Wigner's 
3j coefficients. Our aim below will be to show that the contribution of each 
of these components is determined by its degree of connectivity. Indeed, the 
leading term will be provided by paired graphs 7 € Tp (T) , where the nodes are 
partitioned into disjoint pairs. Next to that, we shall show that the components 

— E 

where at least p nodes are connected are bounded by 0(1^ ) (a bound that can 
be improved for some values of p). This bound can be obtained by partitioning 
these connected graphs into trees, that is subgraphs with no loops, and then 
associating these trees to the coefficients of some unitary matrices arising in 
tensor spaces generated by spherical harmonics. The proof of (|10|l can then be 
simply concluded by a direct graph-counting argument. Let us now make this 
argument rigorous. 

We start by introducing some notation, which is to some extent the same 
as in Marinucci (2005). We need first to introduce a new set of triples C = 



< 



K, 



Pi— pi 



2/n 



(10) 
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{(41,42,4s) , ... (4ii,4i2, 4*s)} , defined by 

(4i, 4-2, 4-3) : =(/n,/i2,ii3) forr=l, ...,2pi 
(4i,42,43) : = (huhziha) for r = p x + 1, ...,2p 2 

(4i,42,43) : = (iji, Il2,li3) for r =pi-% + 1, ...,2p/ ; 

more explicitly, the set £ is obtained by replicating 2pj times each of the 
(ln,li2,ha) triples. Let T be a set of indexes {(r, fc)}, where fc = 1,2,3 and 
r = 1, 2, Y^i=i Pii-f° r anv 7 ^ T(T), we can define 

<5( 7 ;£):= [] (-1)-—^'"" ! (") 

for brevity's sake, we write £(7) rather than £(7; C) whenever this causes no 
ambiguity. Recall that 



Ea 



Z ri k 1 m rikl a l r2k2 m r2k2 — ( 1) rihl Ct kl &i^J>m T ^ h \ 2 



(12) 



In view of ©, and because the spherical harmonic coefficients are (complex) 
Gaussian distributed, the following formula holds: 



A n 



(r,k)£T 



7er(T) 



(13) 



Write {(r, A;),.} € 7 to signify that the pair {(r,k),(r r ,k')} belongs to 7, for 
some (r', fc'); for any diagram 7, we can hence define 



D( 7 ) 



We define also 



e n 

m rh =-e rk r:(r,fe)eT 
{(r,fe),.}£7 



4i 42 43 

m r \ m r 2 rn r 3 



*(7) 



(14) 



0(A) := £i>( 7 ) = £ e n 



7€A 



7eAm ri = -* ri r:(r,k)GT 
{(r,fc),.}ey 



4 4 4 

m r i m r 2 m. r 3 



-5(7) ; 



in words, -D(.) represents the component of the expected value that corresponds 
to a particular set of diagrams. Notice that 



HikU= e e n 



m rfc = -£ rfc I r:(r,fe)£T 
r:(r,k)eT 



4l 42 



tj-3 



m r i m r2 m r3 /11 /c„ 

fc=l v 1 



n 



e e n 

•er(T) m rfc = 
£>[r(T);£] 



7 er(T) m^, r:(r,fc)eT 
{(r,fc),.}e-y 



^rl ^r2 ^r3 

m r i m r 2 m r 3 
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Now 



D[T(T);jC}=D[T p (T);jC}+D[T(T)\T p (T);C} ; 



by an identical combinatorial argument as in Marinucci (2005), it is simple to 
show that 

D[r p (T);£]=n{(2 P< -l)!!A« al8 } . 
To complete the proof, it is then sufficient to establish that 

D\r(T)\r P (T);£] = o((2i 11 + i)- 1 ) . 

It is shown in Lemmas 3.1-3.3 in Marinucci (2005) that diagrams with a 1- 
loop (7 G r^(!)(T)) correspond to summands identically equal to zero, whereas 
diagrams with p nodes and loops of orders 2 or 3 can be reduced to terms 
corresponding to diagrams with p — 2 nodes times a factor 0((2lu + l) -1 )- In 
the sequel, it is hence sufficient to focus only on graphs which have no loops of 
orders 1,2 or 3. 

Now call R the set of nodes of the graphs, and partition it into subsets such 
that 

R = Ri U R 2 U ... U R g . 

Then it will also possible to partition 7 into subdiagrams 7l , 7 2 , j g , 7 12 , l g -i, g 
such that -f 1 includes the pairs with both row indexes in R\, -f 2 includes the 
pairs with both rows in R 2 , 7 12 includes the pairs with one row in R\ and the 
other in R 2 , and so on; we assume all internal subdiagrams 7 f to be non-empty, 
whereas this need not be the case for 7i -. In terms of edges, j x includes the 
edges that are internal to R\, -f 2 includes the edges that are internal to R 2 , 7 12 
includes the edges that connect R\ to R 2 , and so forth. Note that, 



7 = ( 7l U 72 U ... U 7 3 U 7l2 U ... U 7g „ lj9 ) 



and we can write 



9 9 



d m - EE E n 

i=l j=i+l m Tk = ~l Tk I 1=1 
{(r,fe),.}e 7i . 



y TT l ri t r2 c r3 - s{ 

J- J- y m r i m r2 m r 3 

m r . k = -e rk reRi 
.{(r,fe),.}e 7i 



9 9 z Tk f g 1 

EE E n**™ 

i=l j =i+l m Tk = -l Tk \i=l ) 
{(r,fc),.}€7« 



where 



Xn^--= E n 



m rk = —l rk rGRi 
{(r,fc),.}G 7l 



W2 



tr3 



m r i Tn r2 m r 3 
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Xfl i;7 . can be viewed as a vector whose elements are indexed by m Ti ^ i , where 
ri e Ri and {(r^, k,j), .} ^ 7 i (indeed those indexes m ri fc i such that {(rj, fcj), .} e 
7 i have been summed up internally). For instance, for g = 2 we have 



D( 7 ) = 



E 



m r f t ——£ r j t 
{(r,fc),.}G 7l2 



e n 

m r k=—trk reRi 
{{(r,k),.}ey 1 



x < 



^rl ir2 £r3 | x 

m r i m r2 m r3 ' v 11 



y TT ^ ^ <r3 | %) 

J- J- V m r i m r 2 m r 3 

m rk =—£ r k rER2 
{{(r,k),.}e-y 2 



> *(7i2) ■ 



In Figure II, we provide a graph with eight nodes = 8 (right), and then 

(left) we partition it with g = 2, #(i?i) = #(i? 2 ) = 4; the nodes in R\ arc 
labelled with a circle, the nodes in R2 are labelled with a cross, the edges in ■y 1 
and 7 2 have a solid line while those in -f 12 are dashed. Here we have 3 + 3 = 6 
internal sums and six external ones. 



Figure II 

Assume now that -f i docs not include any loop, for i = 1, ...,g. Our point 
will be to show that 

1^(7)1 <nii x ^ 11 as) 

where |.| denotes Euclidean norm, and 

11^7*11 < IT (2£,. fe + l)- 1/2 <(2 min l rh + l)(#(«<)-i)/2 (16) 

{(r,fe),.}€7i 

note that if 7i does not include any loop the number of edges it contains must 
be identically equal to #(i?i) — 1, where #(.) denotes the cardinality of a set. 
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Let us consider 11150 first. It is clear that we can choose new indexes such that 



Xr 4;7 . =: I™ is a vector with elements 



mil ,■ ■ ■ ,771*' 



}■ 



,3 ; 



we write T for this new set of indexes, namely 

(1,1) (Mi) 

(5,1) (9,Vg) 



T 



here g can be viewed as the number of trees and Uj as the number of vertexes 
which are in a given tree i. Clearly the vectors ...,X^ have dimensions 

= IljLi( 2 % + !)■ The following Lemma can be viewed as an extension 
of the Cauchy-Schwartz inequality. 

Lemma 3.1 (Generalized Cauchy-Schwartz inequality) Let 7 be a partition of 
T with no flat edges. 



nn e nk 



mil ■ 



\&m <n|k w 



i— 1 j — 1 m,ij——l 



and 



nn E 



E ••■ E 



i=l j=l m; 



Proof The result follows from the iterated application of the Cauchy-Schwartz 
inequality; we shall argue by induction. Without loss of generality, we can 
assume that the diagram 7 is connected (if it is not, argue separately for the 
connected components). It is trivial to show that the result holds for g = 2, 
indeed in that case it just the standard Cauchy-Schwartz result. Let us now 
show that if the result holds for the product of g — 1 > 2 components, it must 
hold for g components as well. Recall we consider diagrams with no fiat edges, 
so the indexes cannot match on the same vector X^ l K Relabel terms so that 
there exist (at least) a link (that is, a common index m) between the first two 
vectors X^\X^ 2 '. Without loss of generality we can order terms in such a 
way that the matching is internal for the first v* indexes and external (that 
is, with the remaining nodes (3,4, ...,<?)) for my : j = v* + 1, ...,v\ and mij ■ 
j = v* + 1, ,..,V2- We write also 7 = 7 12 U 7^, where 7 12 is the set of edges 
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linking node 1 to node 2 and 7^ = 7\7i2- We have 



9 vt Ui 



g 



nn e \ \u\ x ^,..m iVi \m\ 

i=l j=l mij=—iij 



< 



E ■•• E 



r i 

E "' E [ Y^m\ 1 ...m 1 v 1 ^m\ 1 ...m2V 2 \ I ^(Tl2 ) I 



9 Vi tij 



nn e \\u M ff ^i 

i=3 j'=l mij = -lij 



J=3 



< 



E 

1— 

E 



E 
E 



E ■■• 

mil— —^11 WHv* — — ^1 

1/2 



/ y f | yl mn...ini»i 



1/2 



m2i...m2i)2 



m 2 i = -«2l m 2t ,*=-<;2^* 



9 Vi 



nn e ni^, 

i=3 J = l m,ij=— lij J Ci=3 



1^(712)1 



(17) 



The last step follows again by standard Cauchy-Schwartz inequality. Now 
define 

1/2 



m lv * +1 ...m 2 V2 



mil — —in rn iv* — ~ *--iv* 



E 

Ll — _ 

E 



V \ IxV 



/ y f \ TO21 ...m2«2 



so that H17[l becomes 



m.21 — — £21 m 2v 



1/2 



E - E E - E ra^^IIK 



ll 1 „* + l = -*!l„* »Tl2u2 =_ * ! 2 l ,2 m 31 = -'31 "' . 

< 



gv g *-gvg 

g , 



i=3 



\s(rhs)\ 



np 



i=3 



(18) 



by the inductive step (indeed within the curly brackets we have the product of 
g — 2 + 1 = g — 1 components) . Now notice that 

■ 1/2 

V V (xW 

1 1,d*+1 = -^1u* Tl2i,2 =_ ^2i,2 
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E 



~^2 til r 

y y ) ^ y ^ y [ ^ ^"1] 1 . . . 771 1 »1 J 



n>i,»' + i--'i"' "72i. 2 



mil — — tn mi t 



*21 

E ■■• 

17l21=-^21 m 2 « 



/ y [ ^ 1 mai...ni2i)2 J 



1/2 



1/2 



/ ^ / ^ | yi mn...miui 



mil— — til TTllu 



<21 

E ■■■ 

17121=— ^21 77121,,=— *2 



1/2 



y [x^ 

/ y ^ 1 m 2 1...7Tl27'2 



and thus by substitution into 118|) the proof is completed. 



□ 



Remark 3.2 We provide an example to make the statement of Lemma 3.1 more 
transparent. Take 



X^ X^ X^ — ]■ < m < / i i — 1 9 3 

A mnmi2 i ^ i m 2 im 2 2"i23 ' ^ 1 m3l"i32"i33 ' 1 — l J — 1 ' 'iJ _ 



here 



(1,1) (1,2) 
T={ (2,1) (2,2) (2,3) 
(3,1) (3,2) (3,3) 

Now take for instance 

7=[{(1,1),(2,1)},{(1,2),(3,1)},{(2,2),(3,2)},{(2,3),(3,3)}] , 
and consider the sum 



<frfty>i* (i) * (2) 



771217712277123 771317713277133 



IWt)I 



I /— 1 j—l m 

/ y / y / y / j I 777.1177112 777-1 
77111 77112 "122 771 2 3 

< llx^ll Ik^ 



X (3) I 

171117712277123 771127712277123 I 



by Lemma 3.1, where 



^ ' (-^7711177112) 



E (* 



(2) ^2 

7Ti2l7Ti22"l23 I 1 



7TI21ITI22ITI23 



X (3) 



E (* 



(2) 



m3im 3 2m33 j 



m3im32m33 
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It is not difficult to see that i|15fl is an immediate consequence of Lemma 3.1; 
in particular, note that 7 can be viewed as the diagram which is obtained from 
7 by identifying all nodes that belong to the same set Ri, for i = X,..,,g. Now 
let us consider ifTBj l. The following result uses the previous inequality to bound 
the components of a connected graph. Without loss of generality, we re-order 
terms so that for all r, we have £ r i < i r 2 < i r 3 and 

in < in < ... < im ■ 

Note that the same inequality need not be satisfied for the sequences {i r k] r=1 ^ 
k = 2. 3. 

Lemma 3.2 

a) Every connected graph 7 6 T with no loops of orders 1,2 or 3 can be 
partitioned as 

7 = (71 U 72 U ... U 7s U 7l2 U ... U 7fl _ 1)S ) 

where r ) i has no loops of any order and is such that > 2 (in other words, 

7 can be broken into into binary trees with at least two nodes) 

b) For all j i} i = l,...,g we have 



11**^11 < TT (24 fc + l)- 1/2 <(2 min £ rk + l)(#<*M>/a 

{(r,fc),.}67i 

that is, every tree with p nodes corresponds to summands which are 0((2in + 
l) - ( p-1 )' 2 ), where p is the number of nodes in the trees, 
c) For all connected graphs 7 with p nodes, we have 

p/4 

iDlrXW^HWri + ir 1 . (19) 

r=l 

Proof 

a) We drop edges till we reach the point where there are only binary trees 
or isolated points. Each of these isolated points can be connected to either 
another isolated point, in which case we simply have a tree with two nodes, or 
to a binary tree. The graph which is obtained by linking this point to the tree is 
itself a tree if its paths have at most two edges: recall there are no loops of order 
2 or 3. On the other hand, if the graph has a path which covers four nodes, 
then we delete one edge and obtain two trees with two nodes. These procedures 
can be iterated until no isolated point remains. 

b) We recall the identities (see VMK, chapter 8) 

il i-2 i3 \ ( iA+m 3 1 fji3,ma 

mi m 2 -m 3 ) ' ^2i 3 + 1 ll ~ miA - m2 

r la,m a _ 1 i^-ij+ras./o* 1 1 / ^ ^ ^ \ 



mi mi — ^3 
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where the coefficients C^'^ # 2 m2 (Clebsch-Gordan coefficients) are the ele- 
ments of a unitary matrix which implements the change of basis from a tensor 
product to a direct sum representation for a space spanned by spherical har- 
monics. More precisely, denote by {Yg x ® Y£ 2 } mi m2 the elements of a basis for 
the tensor product <g> Y^ 2 ; here, Yt is the vector space generated by Ye m , 
m = —£, ...,£: see for instance Vilenkin and Klimyk (1991), Chapter 8 for a full 
discussion of tensor products and their properties. For our purposes, it suffices 
to recall the identity 

£ 2 +£i t 

l=\l 2 -l x \m=-t 

More compactly, we might refer to the (2^ + 1)(2£ 2 + 1) x {2£ x + 1){2£ 2 + 1) 
matrix C, whose elements |c'^ T ^ li ^ 2 TO2 | arc indexed by mi,m2 over the rows 
and m over the columns. The matrix C is unitary: it transforms a basis the 
orthonormality relationships read 

E(~<t,m s-il' ,m' X fJ X ml 

L '/i,mi,<j,mj t; /i l mi,<i,mj ~ °m > 

mi,m2 

Ef~tl,m /-tt,m X m iX m ' 2 

U <i,mi,< 2 ,m 2 U <i,Tn' 1 ,<2,ro2 ~ m i m 2 ' 

Now the argument can be iterated to higher-order tensor products, to obtain 

£ 2 +£i i 

<=|< 2 -<i|m=-< 

^ ^ ^ ^ ^J,roi,<2,mj^,m,4,m 3 ^'' 

€=|< 2 -<iK'=K3-<| m=-<m'=-<' 

The orthonormality conditions now read 

f £ , , I 2 f * , , I 2 

El \ " f~,l,m fit' ,m' I _ \ " I \ " s-,£,m fit' ,m' I _ -, 

\ 2-^i U <i,mi,<2,roj u <,ra,< 3 ,m 3 f 2^ | L ^i,roi^ 2 ,m 2 L ^,m,£ 3 ,TO 3 f 1 ■ 

mim2iH3 lm=— J Wm 1 lm= — £ ) 

More precisely, the coefficients 



s^£,£' ,m V f^£,m s~i£' ,m' 

L/ £ 1 ,m 1 ,£ 2 ,m 2 ,£ 3 m 3 '~ 2—1 £l,mi,£ 2 ,m 2 { ~'£,m,£ 



^3, "*3 

are the elements of a unitary matrix with (2£\ + 1)(24 + 1)(24 + 1) rows and 
columns; the columns are indexed by £,£',m and they are 

£ 2 +£i £ 3 +£ £ 2 +£i 

]T J2 (2^'+l) = (24 + 1) (2*+l) = (2*i + 1)(24j + 1)(24 + 1) • 

t=\t2-ti\t'=\t3-t\ £=\£ 2 -£i\ 



16 



In general, we have the coefficients 



a 



Ai , A2 , . . . , Ap_i 
£1 ,mi ;.. .;£ p m p 



such that 



Mi= — Ai (j, p2 = — Ap_ 2 



Ap-i 



x...xC* 



A p -i,/j 

Ap-2,Pp_2^P, m P 



E\ "* \ \ "* ^Ai^a.-.^Ap-i^^ 

••• Z^ ^i,mi;...;£pm p J A y _i,/i 

Ai Ap_2 Ap_i /^— — Ap_i 



g> ... <g> y>„) 

l 11 *p J m,...jn 

whence the ortho-normality conditions yield 

2 Ap_i 

Ef/iAi,A a Ap-i;,ul _ r r ,Ai,A 2 ,...,Ap_i;/j , | _ 1 

\ Lj £i,rn 1 ;...-e p m p J Z^"'Z^ Z^ Z^ \ ii ,mi ;.. . \l p m p J ~ 1 ■ 
mx,...m p Ai Ap_2 Ap-i ^i— — Ap-i 

It follows that any cross-product of p Wigner's 3j coefficients corresponding to 
a binary tree can be bounded by a term like 



Ai Ap_i 

E E - E 

m 1 ,...m p ;m fj, t =— Xi V p -i= — A p -i 
•\ -1 

p-1 

(2£+l)Y[(2X J + l 



£1 ^2 Ai 

mi TO2 /ij 



Ai £3 A 2 

Ml "13 A*2 



E 



Ap-i,/Jp_i;^ 



Ai Ap— 1 

E ■■ 

Mi — — Ai Mp-i — — Ap-i 

\ -1 / 

(m + 1) n (2A, + 1) I E «£;:::fc Am } a - 1 IK 2 ^ + 



mi,. . .m p ;m 

p-i 



mi 7 ...m p-.m 



3=1 



c) Connected components with loops of order 1, 2 or 3 can be reduced as 
shown in Marinucci (2005). From part b), we have 



(p-ff)/2 

\D[ T ,C}\< J] £i rl + i)-\ 



(20) 



whereas from a) we learn that it is always possible to choose a partition such 
that g < p/2; the result follows immediately. 

□ 

Remark 3.3 An inspection of our argument (especially in c)) reveals that the 
bound Ijl9(l can be improved in those cases where it is possible to partition the 
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graph 7 into a minimal number of trees. For instance, it is known (see for 
instance Biederharn and Louck (1981)) that all connected graphs with up to 8 
nodes can be partitioned into two binary trees; for such cases, we hence obtain 
the bound 

p/2-1 

\D[ r ,£]\< n (2^i + ir\ P<8. (21) 

r=l 

For p > 10 it is no longer true that such a partition necessarily exists; we leave 
for future research, however, to ascertain whether (|20|l can be improved by a 
more efficient use of graphical arguments. 



We consider now the case where the angular power spectrum is unknown 
and estimated from the data. Define 

uim ■= —f= , uim := —j= , m = 0, 1, I ; (22) 
\ C 

from Marinucci (2005) we recall the following simple result. 



Lemma 3.3 Let I and p be positive integers, and define 



fe=i 

For u and u defined by 11' 21) . we have 



E { Uiq. . . uipui i . . . Uj iu* n . . .u*j . . . ui k . . . Uj k u* k . . . 

qo times (ji times times Qk times q ' k times 

E { u m .„ui un . . .UiiUg . . .U^ .- Uik . . Mi k u* lk . . m^ > x g(l; q + qi + ... + q' k ) . 

q times qi times g J times 9fe times q ' k times J 

Now partition the indexes {£u , £12, £13 , — , £ri , £r2 , £r3 } into equivalence classes 
where l r k takes the same value; we can then relabel them cts £-u , u — 1 , . . . , c where 
c is the number of such classes, i.e. the number of different values of £; clearly 
c < 3i?. For convenience, we write 

[{O11.i12.i13)! -0n ) ii2,i/3);Pi, -,Pi}} ■= {^11,^12,^13, ■■■,£ri,£r.2,£r3} , 
and define 

G[{0ii,ii2,ii3),...0n,ix 2) i/3);pi,...,M] : =f[.9(^;^) 

u=l 

c (#H)/2 ' 2L + l 1 



n n i2~ e +2j _ 1 ) > 

u=l 3 = 1 v " J J 
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where denotes the cardinality of the class where £ u belongs, or in other 

words the number of times that each particular value is repeated in {£n, £12, £13, --,£ri,£r2, £r3}- 
We give some example to make the previous definition more transparent. Con- 
sider for notational simplicity the univariate case 1 = 1. For h < h < h we 
have c = 3 and 

3 R 

g [{(in, in, 113), --(111,112, ii3)-,pi, -;vi}] = n n 

u=l j = l 

for h = I2 < h we have c = 2 and 



2/ u 



91 



2 j - 1 



G[{(in,W U )^ ; 
for li = I2 = ^3 we have c = 1 and 

3i? 

G [{(I11J12, 113), --(111,112, li3);Pi, --,Pi}} = Y[ 

i=i 

It is shown in Marinucci (2005) that (see J22J0 

E^ l2l3 =EI% 2h g(l u l 2 ,W,p)-, 

for any peN; likewise, take pi, .-.,pi £ N: by exactly the same argument it is 
immediate to obtain 

(24) 

Notice that. 

lim G[{(ln,li2,hs),— (hi, I12, ha); Pi,— ,Pi}] = 1 • 

Zll — s-OO 



2*1 + 1 
2/i + 2j - 1 



The following result extends Theorem 3.1 to the case where the angular power 
spectrum is estimated from the data. 

Theorem 3.2 Assume that (^1,^2,^3) 7^ (h'l, h'2, h'3) whenever i ^ i' . There 
exist an absolute constant K Pl _,_ PI such that, for pi > 1 , i = 1, 



i=l 



< 



K 

Ha + 1 



Proof It is sufficient to combine Theorem 3.1 and (|24|l . 



□ 
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An immediate consequence of the Theorems 3.1 and 3.2 is the following 



Theorem 3.3 (Multivariate Central Limit Theorem) For any k e N, as In — > 

oo, 

( j[lljl2jl3 Ilkllk2lk3 J | ^11^12^13 ^fclifcjjfcg | ^ 7V(0 Jfc ) 

\f^hlll2ll3 ' ' V ^ l lllk2lk3 J ' V V Ai n'l2'l3 ' ' V^hllk2lk3 J 

where denotes the (fc x 1) identity matrix. 

Proof The results follow immediately from Theorems 3.1, 3.2 and the method 
of moments. 

□ 

Remark 3.4 It is immediate to see that, for some K' ^K'^^ > 

KL 



\G[{(hl,h2,ll3),-(lll,ll2,ll3);Pl,-,Pl}] - 1| < 



pi---pi 
2/ii + 1 

whence Theorem 3.1 can be also formulated as 



i=l 



if" 
" 2/n + 1 



A careful inspection of the proofs reveals that the rates of the bounds are exact 
(we shall come back to this point in the next Section); in other words, the 
bispectrum converges to a Gaussian distribution with the same rate when either 
the angular power spectrum is known or unknown (the bounding constants 
differ, however). This result settles some questions raised in Komatsu et al. 
(2002)), where the distributions of Ii n i i2 i i3 and Ii n i i2 i i3 where compared by 
means of Monte Carlo simulations. 

Remark 3.5 It is interesting to note that the angular bispectrum ordinates at 
different multipoles are asymptotically independent, for any triples (In, I12 , Z13) 7^ 
(Z21, I22, ^23)- 



4. HIGHER ORDER RESULTS 



The results of the previous sections can be further developed to provide 
some higher order approximation on the moments of the angular bispectrum. 
For notational simplicity, in this Section we shall focus on the univariate case, 
that is, we do not consider cross- moments; the multivariate generalization is 
straightforward, however, and no new ideas are required. In the statement of 
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the Theorem to follow, we use the Wigner's 6j symbols, which are defined by 



a b e 
c d f 



a,/3,7 e : S,cj> 



a b e \ I c d e \ f a d f 
a (3 e J I 7 5 —e ) \ a 5 —(f) 



(25) 



see VMK, chapter 9 for a full set of properties; we simply recall here that 

1 1 1 



a b e 
c d f 



< min 



V '(2a + l)(2c+l)' v/(26+l)(2d+l)' y/(2e + l){2f + 1) 

(26) 



Theorem 4.1 There exist an absolute constant K p such that 



EIl " hh - (2p- 1)!! - f^-^Wi.Ms) !> (2y<- I >" 



where 



G 



6 6 6 f l! / 2 Z 3 

2ii + l 2l a + l 2fe + l 1 *i *2 ? 3 



< 



(2Zi + l) 2 



(27) 



, for h < h < h , 



96 



24 



2l 2 + 1 2Z 3 + 1 
96 24 



22 2 + 1 2/i + 1 



h 


h 


h 


h 


h 


h 


h 


h 


h 


h 


h 


h 



, for l x = l 2 < h 
, for h <l 2 = h 



ka{1, I, I) := 



6 x 18 2 



21 + 1 



I I I 
I I I 



, for h = h = h = I 



It holds that (see 



C 
2h + 1 



some C > 



Proof For all p <E N, a direct combinatorial argument yields the following 

^W aIs = (2p-l)»Af ll3la (=:fli(ii,Z 2 ,M) 
+ ( f W^a.Wp-BjHAfJ*,} (=:D 2 (Zi,i a ,i a )) 



+ - 



2p \ ( 2p-4 



2 l 4 



Z 2 ,/ 3 )(2p-9)!!Af i - 4 Z3 



2p 
6 



K6 (/ 1 ,Z 2 ,; 3 )(2p-7)!!Af i - 3 i3 



(=: D 3 (h,l 2 ,l 3 )) 



21 



+ 



+i 

2p 
6 



?)(? 

2p- 6 
4 

+ 



4 



2p- 8 
4 



^(/l.fc.is) x (2p-13)!!Af i - 



-6 

■2h 



2p- 10 
2 

2p 



K6(/l,i2,/3)K4(^l,i2,/3) X (2p - 



11)!!A?J 2 5 ;3 



K8(/i,i2,/ 3 )(2p-9)!!Af i 7 



+... + K2 P (li,h,h) {=■ D p (l u l 2 ,l 3 )) , 

where the sum runs over all elements such that the factorials are nonnegative; 
also, we take (2p — k)\\ = 1 when k > 2p and 

K u (h,l2,h) ■= J~]D[T,h,h,h] ; 



► (=:D4(Jl,i2,i3)) 



where the sum runs over all possible connected graphs with u nodes, that is 
Ku{h,h,h) represents the expected value corresponding to connected compo- 
nents with u nodes. The term Di(li,l2,h) corresponds to the sums over all 
graphs with exactly p connected components, each of them with exactly two 
nodes. The term -D2C1, ^2, ^3) correspond to the sums over all graphs with ex- 
actly p — 1 connected components: the number of such graphs corresponds to 
the possible ways to select 4 nodes out of p and then partition the remaining 
2p — 4 into pairs, that is. 



2p 
4 



(2p-5)!!Af; 



2 

hh 



The argument for the remaining terms is entirely analogous. Now from the 



min(2,f). 



for u > 6, whence 



previous section we know that K u (l\, I2, h) = 0(l 1 
the result will follow from an explicit evaluation of K4(Zi, h, ?3)..For the latter, 
note that connected graphs with four nodes must include a loop of order 1 , 2 
or 3. Terms with a 1-loop are identically zero (Marinucci (2005), Lemma 3.1). 
Graphs associated to terms with 2- or 3-loops (7 G ^L(2) or 7 e ^l(3)) are 
represented in Figure III. 



Figure III: 7 e r L(2) ,r L(3) 
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We consider the three cases separately. 

a) For li < l 2 < h, we consider first the graphs 7 G ^L(2) (on the left- 
hand side of the figure). By means of Lemma 3.2 in Marinucci (2005), it is 
simple to ascertain that each corresponding term produces a factor (2Zj + 1), 
where the k correspond to the index which is not in the two-loop (the dashed 
line) . Direct counting of permutations shows that there are 6 such graphs for 
each fixed £j. Likewise, the Wigner's 6j coefficients arise in connection with the 
graphs 7 G r L ( 3 ) where each node is linked to all three others, compare l|2"5|) 
and Lemma 3.3 in Marinucci (2005); direct counting of possible permutations 
shows that there can be six distinct combinations of this form (fix for instance 
node 1, which by assumption is linked to all the other three nodes: there are 
three degrees of freedom to choose the connection with node 2, then two left for 
node 3, for a total of six, as claimed). 

b) For l\ = I2 7^ £3, there are two possible types of graphs 7 € ^L(2), that 
is, those where both 2-loops involve I2 and those where one 2-loop involves h 
and the other I3. In the former case, there are three ways to choose the pairs, 
two ways in each pair to choose the links, and two ways to choose the way to 
match the edges corresponding to I3: the total is 24. In the latter case, it can 
be checked that there are 6 possible choices of pairs and 4 possible matchings 
for the 2-loop which involves I2 and Z 3 . The remaining term is similar. 

c) For lx = l 2 = £3, again we choose in three possible ways the pairs, plus the 
single link in two possible ways; within each pair, we can choose 3x3 couples 
and two possible ways to link them. The remaining term is similar. 

□ 



Remark 4.1 Theorem 4.1 can be used to establish (in a merely formal sense) 
Edgeworth or Cornish-Fisher type approximations for the asymptotic behaviour 
of the angular bispectrum (see Hall (1991)). We provide here only some heuristic 
discussion, and leave for future research the possibility to establish rigorously a 
valid Edgeworth expansion. Let 



/ 7-2P 
(2P) ._ E I I hhl 



'3 



A? 



hh. 



whence 



£<U^ - E^^-^ + ^^E^(2p-2)(2,-l)!! + 0(Z-) 

p=l P ' p=0 y ' p=2 F ' 



p=0 P' p=0 p ' 



em— 
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Again in a formal sense, we can take Fourier transforms on both sides, leading 
to the conjecture that 



Pr 



h 



I 1 ' 2^3 



< x > = 



(2/3 



(p(z)dz- 



24 



(z i -6z+3)<f>(z)dz+0(l^ 



where <j){z) denotes the density function of a standard Gaussian variable. It 
is remarkable that the first term in this expansion involves the fourth Hermite 
polynomial H±{z) — (z 4 — 6z + 3) rather than the second Hi(z) = (z 2 — 1) as it 
is more commonly the case. 

Remark 4.2 By means of (|24|) . it is immediate to extend Theorem 4.1 to cover 
the case where the angular power spectrum is estimated from the data; we 
provide below a detailed calculation for 2p — 4. Also, higher order moments can 
also be evaluated iteratively, according to the following expression which holds 
for p > 3 

EI% 2l3 =(2p-l)llAl l2h 
f )(2p-5)!!Af- 2 i3 x[^ 2i3 -3A^ 3 ] 

f ) (2p- 7)!!Afj 2 3 i3 x [EI? ll2l3 - 15Af lW3 ] + ...+ 



2p 

2p-2 
(2p-l)!!Af iW3 +£ 

fc=2 



^/Lt ~ (2? - 3)!!Af^ 3 l + K2p(ii,i 2> l 3 ) 



2p 
2fc 



(2p-2fc-l)!!AfJ* 8 X [25J$ aI , - (2fc - l)!!Af i/2i3 ]+ K2p (M 2 ,*3) 



This expression can also be exploited to derive approximations of moments, 
where the term K2 P (h, h, h) is simply neglected. 

Remark 4.3 It is interesting to note that for p = 2 H27JI holds with K p = 0. We 
provide hence an explicit evaluation of these moments. First recall that (VMK, 
Chapter 9) 



^1 h h 
h h h 



{(h + h~ h)Kh + h- I2W2 + h- hV-Y 



-1 2 



(h + h + h + 1)! 



2l 2 +2l 3 

E 

n=h+l2+h 



(-l) n (n+l)! 



((n -h-h- hY-) 4 mi + 2l 2 - n)\)((2h + 2l 3 - n)!)((2I 2 + 2l 3 - n)!) 
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This expression can be simplified for some values of the triple (li,h,h)- More 
precisely, we have 

I I I I _ A (-!)"(» + !)! 
* * 1 J ((n-3*)!) 4 ((4Z-n)!) 3 

and (VMK, eq. 9.5.2.10) 

Zi Z 2 Zi+Z 2 1 (2*i)!(2* 2 )! 



*i *2 Z1+/2 J (2Zi + 2Z 2 + 1)! 



Thus we have 

"' + 2*i + l ' 2/ 2 + l ' 2*i + 2* 2 + l ' "(2Zi+ 2Z 2 + 1) 



FT* -3 + 6 + 6 6 , fi (Wj , 



(28) 



,4 _ Q .o2 , 96 , 24 , _ ((2*)!) 2 



£1,1 2 , = 3 x 2 2 + + + 48 ;\ J '\, , (29) 

u ' 21 21 + 1 4* + l (4Z + 1)! y ' 

and 

F r4 _ , v fi 2 + 6X182 , fi 4 a') 6 ^ (-!)>+!)! 

The validity of l|28|) - (|30[l has been confirmed with a remarkable accuracy by 
a Monte Carlo experiment (not reported here), where 200 replications of the 
sample bispectrum at various multipoles (h,h,h) were generated for Gaussian 
fields; their sample moments were then evaluated and found to be in excellent 
agreement with our theoretical results. 

By means of (|24[l these results can be immediately extended to the case 
where the normalization is random. More precisely, we obtain for h =/= I2 



6 6 6 (2*01(212)1 

V JT 2/1 + I 2/2 + I 2/1 + 2/2 + I (2/i + 2*2 + l)! 
2*i + 1 2*2 + 1 2*3 + 1 1 



2*i + 3 2*2 + 3 2/3 + 3 _ 



EI* = 02 + xr^ + -^+48 



96 24 ((2*)!) 2 



and 
El, 



in 



2Z + 1 4Z + 1 (4/ + 1)! 
2* + 1 2* + 1 2* + 1 4* + 1 
2Z + 3 21 + 5 21 + 7 4* + 3 



o ::36 6xl8 2 (*!)6 * (-!)"(» + 1)1 

2* + 1 + ((3* + l)!) 2 ^ ((n - 3Z)!) 4 ((4* - n)!) 3 

,; ' 2* + l 



n 



X 

fc=l 



2* + 2fc - 1 
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These expressions can be of practical interest for statistical inference on cos- 
mological data. For instance, the square bispectrum is often used in goodness- 
of-fit statistics to test the validity of the Gaussian assumption. The previous 
results yield immediately its exact variance, which so far has been typically 
evaluated by Monte Carlo simulations (see for instance Komatsu et al. (2002)). 
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